Scattering tightly bound dimers off a scattering potential 
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The motion of two attractively interacting atoms in an optical lattice is investigated in the presence 
of a scattering potential. The initial wavefunction can be prepared by using tightly bound exact 
two-particle eigenfunction for vanishing scattering potential. This allows to numerically simulate the 
dynamics in the generation of two-particle Schrodinger cat states using a scheme recently proposed 
for scattering of quantum matter wave solitons. 
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I. INTRODUCTION 

For ultra-cold atoms in an optical lattice [EJ3, [H dy- 
namical aspects include transverse resonances [4j density 
waves Q , the evolution of quantum fluctuations @ , the 
speed of sound @, [I[ and time-resolved observation and 
control of superexchange interactions 0] . The aim of the 
present manuscript is to perform exact two-particle dy- 
namics in an optical lattice similar to what has been sug- 
gested in Ref . |ldj , a bright soliton in a one-dimensional 
waveguide. As the dispersion relation for the bound two- 
particle states in the lattice approach case without lattice 
for suitable parameters, this can be used to quantitatively 
test the TV-particle predictions of Ref. [lOj via exact nu- 
merics on the two-particle level for which a soliton is 
simply a dimer. 

Besides the analytical TV-particle quantum mechani- 
cal calculations [10| . the scattering of the soliton has 
also been investigated via numerical methods on the 
TV-particle level [TlJ. Different approaches to obtain 
such Schrodinger cat states or related fragmentations 
have been investigated in Refs. [HI, [l3j]. Contrary to 



Schrodinger cat states of a single atom [14j, cat-like 
states of radiation [l5| or mesoscopic spin-squeezed states 
(which have already been realised experimentally 16]), 
the experimental realisation of Schrodinger cat states 
of say, 100 atoms, is still a challenge of fundamen- 
tal research. Suggestions how interesting quantum su- 
perpositions might be obtained can be found, e.g., in 
Refs. [13, d, Ea HE [SH, H3, [H, H3, HI] and references 
therein. 

For bright quantum matter wave solitons [H, [U HH , 
the mean-field (Gross-Pitaevskii) limit has been shown to 
be achieved already for particle numbers as low as TV gi 
3 [1^]. Many of the pape rs p ublished after the ground- 
breaking experiments [30, [3l| solve the Gross-Pitaevskii 
equation for solitons. However, any mesoscopic entangled 
state which involves superpositions of wavefunctions can- 
not be described by a non-linear equation and therefore 
the reasoning of Ref. [29| is not valid in the situation 



considered here. Thus, instead of applying the Gross- 
Pitaevskii equation, the TV-particle Schrodinger equation 
has to be used to reveal true quantum behaviour of a soli- 
ton created from a Bose-Einstein condensate. Under ex- 
perimentally realistic conditions, the Schrodinger equa- 
tion is given by the analytically solvable Lieb-Liniger(- 
McGuire) model. 

The challenge of the generation of mesoscopic superpo- 
sitions via scattering of solitons is that to add a scattering 
potential removes the separability of the centre-of-mass 
motion and the relative motion; in order to avoid that 
the scattering potential acts like a beam splitter on each 
single atom (rather than the entire soliton), the initial 
state has to be prepared carefully. Mesoscopic entangled 
states with the soliton being in a quantum superposition 
with « 50% probability of moving to the right /left should 
thus be obtainable. The probability to find - in a single 
measurement - (at least) one particle moving to the right 
and at (at least) one particle moving in the other direc- 
tion will be negligible. However, this will not be enough 
to prove that the two parts of the wavefunction really are 
in a quantum superposition - if someone claims that a 
coin is in a quantum superposition of heads and tails, an 
experiment showing only the classical outcomes would 
hardly convince anyone. The experimental verification 
could be delivered via interference experiments [To| . 

Rather than dealing with bright TV-particle quantum 
solitons, this paper treats a simpler but nevertheless in- 
structive case: dimers in an optical lattice. The paper 
is organised as follows: after a short summary of how to 
describe the scattering of bright solitons analytically [l(| 
(Sec. [H]), the two-particle bound states used to describe 
the scattering of the dimer are introduced in Sec. lIIII Sec- 
tion IIVI shows the numeric results in the limit where the 
motion in the optical lattice mimics the motion without 
lattice. 



II. EFFECTIVE POTENTIAL APPROACH 
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The Hamiltonian of the Lieb-Liniger-McGuire [32l . |33| 
model with attractive interaction and an additional 
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scattering-potential V is given by 



1=1 



(1) 



Bright solitons [30] are well described by this model. 
For V = 0, exact eigenfunctions of this Hamiltonian are 
known. Solutions corresponding to TV-particle solitons 
with momentum k read: 
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The corresponding energies are given by 
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is the ground state energy of the system [33(] . As long as 
the kinetic energy is not too large, these states are sep- 
arated from the first excited internal state (which corre- 
sponds to one particle having left the soliton) by a finite 
energy barrier cx g 2 N 2 (see, e.g., Ref. (27|V 

Had the scattering potential been a function of the 
centre of mass of all N particles (V = V(££„ x v )), 
the situation would have been easy as the centre of mass 
and relative coordinates then still separate. However, the 
potential in the Hamiltonian (TT|) is given by 



V = Y,V{x v ) 



(6) 



It would nevertheless be tempting to argue that, given 
the fact that the particles are tightly bound, they be- 
have essentially as a single particle and one could thus 
approximate V (x v ) by V (X) and thus 



V « NV(X) 



where 
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is the centre-of-mass coordinate. However, this approxi- 
mation can give wrong results (as will be shown towards 
the end of this paper) and the mathematically justi- 
fied effective potential approach: 



H, 



eff 
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d\ + v eS (x) 



(9) 



has to be used. The effective potential is given by the 
convolution [HI 



V eS (X) = Jd N xS(x-^Y^LiXu) (io) 

N 

x |C^, fc (ai)| 2 £^(^) . 



This approach is valid for sufficiently well behaved po- 
tentials (like a laser focus) and for solitons which cannot 
break apart for energetic reasons (see the paragraph be- 
low Eq. 



III. TWO-PARTICLE BOUND STATES 

Two-particle bound states in optical lattices are in- 
teresting both experimentally (34j and theoretically [34|, 
HE HE J37l [H, ; recently even three-particle bound 
states [iO] have been investigated. Within a Bose- 
Hubbard Hamiltonian. 



lattice 
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one can use exact eigen-states to do the numerics (the 
restriction to negative pair interactions U is not neces- 
sary, however the idea is to discuss a case close to the 
bright solitons created from attractive Bose-Einstein con- 
densates; the creation/annihilation operators of particles 
at lattice cite j are denoted by ct/c.). Rather than us- 
ing the approach via Green's functions of Ref. (34[, we 
proceeded in Ref. (38| along the lines of Ref. [4l[ to show 
in a straight-forward but somewhat lengthy calculation 
that the dimer wavefunctions in a tight-binding lattice 
with band-width W and lattice spacing b are given by 
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exp [ikb(y + fj,)] 
exp [ikb(u + fi)]/V2 



[I = v 



where 



IP 



16 J 2 cos 2 (fc6) 



1^1 

|4J| cos(fcb) 



(12) 



(13) 



the coordinates v and [i label the lattice points on which 
the particles "sit" and as a basis the Fock basis is used 
(i.e. v — \x refers to the Fock state with two particles 
at the lattice-site v). The interaction energy of two par- 
ticles sitting at the same lattice site is Hk < 0; for the 
wavefunction to be normalisable for fixed centre of mass, 
kb < ir/2 (which implies < 1). To avoid a breaking 
of the dimer, the energy 
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FIG. 1: A single particle is prepared in a Gaussian wavefunc- 
tion which initially moves to the right. After the particle is 
scattered off a smooth potential (dashed line), the final state, 
the wavefunction is partially moving to the right and partially 
moving to the left. The data was obtained by numerically 
solving the time-dependent Schrodinger equation. 



should be lower than two times the minimal possible en- 
ergy of two particles: 



E two < -4J 
As for the single particle energy band, 
E onc — — 2Jcos(fc6), 



(15) 



(16) 



the dispersion relation for two particles approaches the 
free-particle behaviour in the limit bk — > which can 
thus be used to model particles not restricted by optical 
potentials: 
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By choosing the parameters to approximate the disper- 
sion relation for a two-particle "soliton" (cf. Eq. @) 
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N=2 soliton 
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K Z k 



(18) 



the effective potential approach of Ref. [l(| (cf. Eq. (fTU 
will thus be valid for small enough b. 



IV. NUMERIC RESULTS 

The initial centre-of-mass wavefunction in Ref. [1(3] is 
suggested to be a Gaussian (obtained by preparing the 
system in a swallow harmonic oscillator potential). If the 
system can be described within a centre-of-mass approx- 
imation, the expected result will always look similar to 
Fig. [TJ For the initial Gaussian to be realisable experi- 
mentally, the length scale on which the wavcfunctions of 
the relative coordinates decays has to be smaller than the 
width of the initial wavefunction. Given the fact that the 
dimer cannot break into two free particles, this implies 




FIG. 2: Particle density g(x) after a dimer was scattered off 
a step potential (|19|) . The effective potential (|1U|1 in the in- 
set explains why resonances can occur. They lead the dimer 
(which for energetic reasons cannot break into two free par- 
ticles) to split into three parts rather than the expected two. 
The initial wavefunction is given by Eq. (|20p with a/ [3 — 14, 
k /P = 0.75 and x Q p = -50; Vo2m/3 2 /h 2 ~ 0.197. The cal- 
culations were done on a 1601 x 1601 lattice. 



that in the final wavefunction measuring one particle on 
one side of the barrier would also lead to measuring the 
other particle on the same side - thus the final state in- 
deed is a Schrodinger cat state. 

Figure[5]shows the result for a two-particle soliton scat- 
tered off a step potential, 



step 







else 



(19) 



The initial wavefunction is constructed as a superposi- 
tion of the eigen-solutions (p~2|) with < k < (5 which (3 
ensures the validity of Eq. (15]) (cf. Eq. (fl8|l): 



y(xi,x 2 ;t = Q) oc / dfcexp(-a 2 (fc-fc ) 2 /2) (20) 
Jo 

X^2,k [Xl ~ X ,X 2 - Xq) ■ 

If a is not too small this leads to a Gaussian centre-of- 
mass wavefunction initially centred at xq- Surprisingly, 
for the parameters chosen in Fig. [5] the wavefunction 
splits into three rather than the expected two parts. If 
one of the particles was measured, the second would be 
much closer to the location of the first than the width of 
the three moving wave-packets. The strange behaviour 
of Fig. [2] clearly demonstrates that the centre-of-mass ap- 
proximation is not valid here: the effective potential (Tl7J|) 
as displayed in the inset of Fig. [2] at least qualitatively 
explains why resonances can occur which make parts the 
dimer wave-function stay at the potential. 

However, rather than being able to qualitatively un- 
derstand the scattering behaviour, the effective potential 
approach can also be tested quantitatively for realistic 
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scattering potentials. If the scattering potential is re- 
alised via the focus of a laser beam, the potential can be 
approximated by a Gaussian or even the potential 



V c (x) 



Vn 



cosh(x/^) 2 



(21) 



which is analytically solvable on the single particle level. 
Figure [3] a compares the wave-function after scattering 
for the exact potential 



V{xi,x 2 ) = V c (xi) + V c (x 2 ), 
the center of mass approximation 

V(x 1 ,x 2 )^2V c {{x 1 +x 2 )/2) 
and the effective potential [l(| 



(22) 



(23) 



V( Xl ,x 2 ) ~ 2V c ((x 1+ x 2 )/2) + V " {{X1 t X2)/2) , (24) 



where f3 is given by Eq. ([3]). The effective potential 
clearly is a huge improvement over the center-of-mass ap- 
proximation. It still remains to be shown that the final 
state indeed is entangled in the sense that whenever one 
measures one particle on the right (left) side, the other 
particle will be on the same side. To do this, a Measure If 
An Outgoing Wave is Unentangled, 



Pu 



x(*) = 



dx 2 \^(x u x 2 ;t)\ 2 



(25) 



where the wavefunction is taken to be normalised to one. 
Figure [3] b demonstrates that for the wave-function to 
be entangled, the initial velocity cannot be too large. 
Experimental verifications that such states are indeed 
quantum superpositions could be done via interference 
experiments [10| . 



V. CONCLUSION 

To summarise, scattering of bright quantum matter 
wave solitons fig . [IiT | has been investigated for two at- 
tractive atoms in an optical lattice via exact numerics. 
It has been demonstrated that the beyond centre-of-mass 
approximation approach of Rcf. [10| is indeed necessary 
to describe the physical situation at least qualitatively: 
the effective potential derived in Ref. [l(| can explain 
the break-down of the centre-of-mass approximation ob- 
served for specific parameters in the numerics (Fig. [2]). 

For a scattering potential given by a laser focus (the 
potential would then be approximately Gaussian) such 
an effect does not occur. The wavefunction then splits 
into two parts as shown for a single particle in Fig. [1] and 
the effective potential even gives excellent quantitative 
agreement (Fig. [3]) with exact numerics. 
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FIG. 3: a) Modulus squared of the scalar product after scat- 
tering of the solution of the approximate Schrodinger equa- 
tion corresponding to Eq. (|23p with the unapproximated so- 
lution (both calculated on a 801 x 801 lattice). The po- 
tential is given by Eq. (|22|l . From top to bottom: i/3 = 
2.5 , 1.5 , 1. Inset: with the approximation expressed by 
Eq. (|24|l . the agreement is much better. The lines are a 
guide to the eye. b) If the soliton splits into two parts, 
Pmiaou ~ (Eq. (|25p ; plotted as a function of the cen- 
tre of the momentum distribution) indicates entanglement 
(for i/3 = 2.5). The wave-packet was defined to be tx 

j(k*/k o )0 dfcexp ^_ a 2 (A , _ k *Yj 2 ) fa ik ( Xl - X0jX2 - Xo )- V 

was chosen such that Eq. (J23J) predicts 50 : 50 splitting; 
fco//3 = 0.75). 
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